A Two-Field Quintessence Model 
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We study the dynamics of a quintessence model based on two interacting scalar fields. The model 
can account for the (recent) accelerated expansion of the Universe suggested by astronomical obser- 
vations. Acceleration can be permanent or temporary and, for both scenarios, it is possible to obtain 
suitable values for the cosmological parameters while satisfying the nucleosynthesis constraint on 
the quintessence energy density. We argue that the model dynamics can be made consistent with a 
stable zero-energy relaxing supersymmetric vacuum. 



Recent observations of type la supernovae, together 
with Cosmic Microwave Background (CMB) and cluster 
mass distribution data jl]-[| indicate that the Universe 
is flat, in agreement with the inflationary prediction, ac- 
celerating and that the energy density of (baryonic plus 
dark) matter is smaller than the critical density. Thus, 
observations suggest that the dynamics of the Universe 
at present is dominated by a negative pressure compo- 
nent, the main candidates being a cosmological constant 
and a slowly- varying vacuum energy, usually referred to 
as dark energy or quintessence (an evolving vacuum 
energy was discussed somewhat earlier, e.g. ||). The 
main difference between the cosmological constant and 
quintessence scenarios is that, for quintessence, the equa- 
tion of state parameter, wq = p/p, varies with time and 
approaches a present value wq < —0.6, whilst for the 
cosmological constant, it remains fixed at w\ = — 1. 

Several quintessence models have been put foward, 
most of them based on a scalar field which was sub- 
dominant in the early Universe and, more recently, has 
started to dominate the energy density of non-relativistic 
matter. Theoretical suggestions include a scalar field en- 
dowed with exponential [[7-10 or inverse power law po- 
tentials [jy], the string theory dilaton in the context of 
gaugino condensation E^] , an axion field with an almost 
massless quark [[l3| , scalar-tensor theories of gravity [jjdj , 
or one of the fields arising from the compactification pro- 
cess in the multidimensional Einstein- Yang-Mills system 
iPI - Some of these models address the "cosmic coinci- 
dence" problem i.e. the question of explaining why the 
vacuum energy or scalar field dominates the Universe 
only recently. In tracker models, the tracker field rolls 
down a potential according to an attractor-like solution 
to the equations of motion, causing the energy density 
of the quintessence field to track the equation of state of 
the background energy component independently of ini- 



tial conditions |Ll). However, in these models, the overall 
scale of the potential has to be fine-tuned in order for 
the quintessence energy to overtake the matter density 
at present. 

In /c-essence models [Q, as a result of the dynam- 
ics, tracking of the background energy density can only 
occur in the radiation epoch; at the onset of matter- 
domination, the fc-essence field energy density drops 
sharply, increasing again and overtaking the matter en- 
ergy density at roughly the current epoch. At least in the 
original proposal, these features require the introduction 
of a non-linear kinetic energy density functional of the 
scalar field and adjusting it to obtain the desired attrac- 
tor behaviour. 

A commom feature of the proposals mentioned above 
is that the asymptotic accelerating behaviour of the Uni- 
verse is driven by the dynamics of a single field. In this 
work, we shall consider instead a two-field model. Two- 
field quintessence models were previously considered, in 
an attempt to explain how to obtain a small but non- 
vanishing cosm olog ical constant |17| and in the context 
of SUSY QCD |[lq| . Actually, there are several motiva- 
tions for studying potentials with coupled scalar fields. 
In fact, if one envisages to extract a potential suitable 
for describing the Universe dynamics from fundamental 
theories, it is most likely that an ensemble of scalar fields 
(moduli, axions, chiral superfields, etc) will emerge, for 
instance, from the compactification process or from the 
localization of fields in the brane in multi-brane mod- 
els or from mechanisms responsible for the cancellation 
of the cosmological constant (see e.g. and references 
therein). Furthermore, coupled scalar fields are invoked 
for various desirable features they exhibit, as in the so- 
called hybrid inflationary models p0| |2l[ an d in reheating 
models in the presence j22j or absence [|23| of parametric 
resonance. Finally, it has been recently pointed out that 
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an eternally accelerating Universe poses a challenge for 
string theory, at least in its present formulation, since 
asymptotic states are inconsistent with spacetimes that 
exhibit event horizons [M. Moreover, it is argued that 
theories with a stable supersymmetric vacuum cannot 
relax into a zero-energy ground state if the accelerating 
dynamics is guided by a single scalar field f24| . In this 
paper, we present a two-field model whose solutions allow 
for, at least, a partial fixing of these inconsistencies. 

Another interesting feature of our model is that it 
presents two types of solutions, namely one in which the 
Universe accelerates forever and one in which it is possi- 
ble for the Universe to exit from a period of accelerated 
expansion and resume decelerated expansion. The lat- 
ter type of solution is compatible with the conceptual 
framework underlying string theory. 

It is believed that scalar fields with potentials of the 
type 

F(^) = e-^P(^</0 , (1) 

where P((p, ip) contains polynomial as well as interacting 
terms in <p and ip, arise in the low-energy limit of funda- 
mental particle physics theories such as string/M-theory 
p5|, N = 2 Supergravity coupled with matter in higher 
dimensions |2q| a nd phenomenological brane-world con- 
structions |19p7|| . The overall negative exponential term 
in <fi signals that this is a moduli type field which has ac- 
quired an interacting potential with the ip field. A simple 
possibility is 

P{4>, ip) = A + - <M 2 + b (V - Va)) 2 

+C ^-^ a ) 2 +D^-M 2 , (2) 

in units where M = (8-kG)^ 1 ^ 2 = % = c = 1. Notice 
that this potential, for B = C = D = 0, in which case 
only the <p field is present, coincides with the one pro- 
posed in Ref. Q, hereby referred to as the AS model. 
We have considered just tree interacting terms as they 
already capture the main aspects of the coupled dynam- 
ics we are interested in. As we will show, an important 
property of model is that it can lead to an asymptotic 
dynamics where either ip or both fields do not necessar- 
ily settle in their minima at present, which is the key to 
evade some of the conclusions of Refs. [ p3| , concerning 
the stability of the supersymmetric vacua. 

We consider a spatially-flat Friedmann-Robertson- 
Walker (FRW) Universe containing a perfect fluid with 
barotropic equation of state p 7 = (7 — l)/0 7 , where 7 is a 
constant, < 7 < 2 (for radiation 7 = 4/3 and for dust 
7 = 1) and two coupled scalar fields with potential given 
by Eq. (|l|). The evolution equations for a spatially-flat 
FRW model with Hubble parameter H = a/a are 



F=--(p 7 +p 7 + 2 +V> 2 ) , (3) 

p 7 = -2,H{p 1 +p 7 ) , (4) 

= -3Hj> - d^V , (5) 

$ = -3HiP - d^V , (6) 



where d^^V = g^ry, subject to the Friedmann con- 
straint 

tf 2 = ^ 7 + ^ + ^ 2 + u) , (7) 

The total energy density of the homogeneous scalar fields 
is given by p Q = 4> 2 /2 + iP 2 /2 + V((p, ip). 

A necessary and sufficient condition for the Universe 
to accelerate is that the deceleration parameter, q, given 
by 

q = ~ = ^(l+3w Q Q Q +Q r ) , (8) 

where Q r is the fractional radiation energy density, is 
negative. 

We now study how the solutions of the system above 
depend on the parameters of the potential and initial 
conditions. We integrate from N=-30, corresponding to 
the Planck epoch; nucleosynthesis occurs around N=-10, 
radiation to matter transition around N=-4 and N=0 to- 
day. 

There are essentially two realistic types of behaviour, 
illustrated in Figures (Model I) and | (Model II). Model 
I (A = 9.5, A = O.O2,0o = 29,^q = 15, B = 0.002, C = 
6 x 10 -4 , D = 4.5) corresponds to the case where vacuum 
domination, which occurs when Qqo > 1/2, is perma- 
nent and Model II (A = 9.5, ,4 = 0.1, cp = 29, ip = 20, 
B = 0.001, C = 8 x 10" 5 ,L> = 2.8) to the case where 
vacuum domination is transient. Permanent and tran- 
sient vacuum domination have also been found in the 
AS model, in Refs. || and p(i|l , respectively. Of course, 
there remains the (non-realistic) case where accelerated 
expansion never occurs. 

In both models the equation of state has reached 
wq ~ — 1 for the present time, as favored by the avail- 
able data |2^] (and making it hard to distinguish from 
a cosmological constant) but, whereas in Model I ujq re- 
mains negative, in Model II it is in the process of increas- 
ing today towards positive values, then oscillates slightly 
until it reaches its asymptotic value. Similarly, in both 
models, the deceleration parameter is negative today but, 
whereas for Model I q remains negative, in Model II it os- 
cillates and becomes positive before it reaches its asymp- 
totic value. 

Permanent vacuum domination takes place when at 
least the (p field ends up settling at the minimum of the 
potential, thus corresponding to a cosmological constant; 
in Model I, both fields settle at the minimum of the po- 
tential, see Figure ^. Transient vacuum domination oc- 
curs either when the potential has no local minimum or <p 
arrives at the local minimum with enough kinetic energy 
to roll over the barrier and resume descending the poten- 
tial. Notice that the evolution of ip is slight compared 
with (f>, especially for Model II. 

Both models satisfy present bounds on relevant cos- 
mological observables. The tightest bound comes from 
nucleosynthesis, Qq{N ~ — 10) < 0.044, requiring A > 9 
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FIG. 1. The upper panel shows the evolution of Qq (solid), 
il r (dashed) and fi m (dot-dashed) for Model I, corresponding 
to permanent vacuum domination. The lower panel shows the 
evolution of wq (solid) and q (dashed) . 




FIG. 2. The upper panel shows the evolution of Qq (solid), 
£l r (dashed) and Q. m (dot-dashed) for Model II, corresponding 
to temporary vacuum domination. The lower panel shows the 
evolution of wq (solid) and q (dashed) . 



p9| . The bound arising from the most recent CMB data, 
£Iq < 0.39 at last scattering, is less stringent than the 
nucleosynthesis bound. Other bounds we take into ac- 
count are: ft m = 0.3 ±0.1 w Q < -0.6, ft Q = 0.65 ±0.05 
|(| and h = 0.65 ± 0.05 § today. Indeed, Model I has 
h = 0.6, tt Q = 0.7, fl m = 0.3, wq = -1 and q = -0.5 
today, £Iq = 0.042 at nucleosynthesis. Similar values are 
found for Model II, namely h = 0.6, £Iq = 0.7, fl m = 0.3, 
Wq = —0.9 and q = —0.4 today, Qq = 0.042 at nucle- 
osynthesis. 

Our models seem to be more sensitive to changes in 
the initial conditions than models with just one scalar 
field and, in particular, the AS model (this is to be ex- 
pected since there is more freedom e.g. in the way ki- 
netic energy is shared between the two fields) but no fine 
tuning of the initial conditions is needed. Indeed, fixing 
e.g. x = z, corresponding to equipartition of kinetic en- 
ergy between </> and ip, we have studied the behaviour 
of pq and seen that, for a wide range of the remaining 
initial conditions, after some initial transient, each solu- 
tion scales with the dominant matter component until 
Pq begins to dominate. 

We have studied the nature of our solutions for a rather 
broad range of parameters of the potential. We have 
found that it is possible to obtain permanent or tran- 
sient vacuum domination, satisfying present bounds on 
observable cosmological parameters, for various combi- 



nations of the potential parameters. 

A relevant issue of our proposal is that it allows evad- 
ing the conclusions of Refs. j24|, in what concerns the 
stability of a supersymmetric potential. The main argu- 
ment presented in |24[| relies on the fact that, in a su- 
persymmetric theory, one expects that the asymptotic 
behaviour of the superpotential is given by W{4>) — 
Woe^ a ^^ 2 , which, in order to ensure the positivity of 
the 4-dimensional potential V{<j>) = 8\d<pW\ 2 - 12\W 2 \ 
implies that |a| > \^6. However, this value is inconsis- 
tent with the requirement of an accelerated Universe at 
present \a\ = ^/3(1 + ujqo)/2 < 1.5 Q as data suggest 
that ujqo < —0.6 pq] . The situation is different in the 
presence of fields that do not reach their minima asymp- 
totically, as in Model II. Indeed, in this case, the asymp- 
totic behaviour of the superpotential would be better de- 
scribed by the function W((j)) = Woe _Q */ 2 F((/>, i/j), where 
F((j>, tji) is a polynomial in the fields 4> and ip- The positiv- 
ity condition now reads: aZ-Q+AKd^F^ + td^F^/F 2 > 
0. One can then easily see that, by a suitable choice of the 
polynomial F((j),ijj), the positivity condition can be rec- 
onciled with the requirement of successful quintessence. 
Furthermore, since in Model II acceleration is transient 
and occurs only at present, this model is consistent with 
the underlying framework of string theory as is does not 
present cosmological horizons that are associated with 
eternally accelerating universes. 
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FIG. 3. Evolution of the quintessence fields, <f) (solid) and 
ij) (dashed), for Models I and II. 



Solutions corresponding to transient acceleration have 
not been found in previous two-field quintessence mod- 
els. In the model of Ref. [jl8j, where the fields invoked 
are the vacuum expectation values of SUSY QCD chi- 
ral superfields, quintessence energy density grows with 
respect to matter as PqI p m ~ a 3 ' 1+r ^ 2 , where r is the 
ratio between the number of flavours and the number of 
colours. Similarly, in the model of Ref. |r7|| , where the 
potential V(a,$) = e" 4C<T (A+ ±m 2 $ 2 [l + 7 sin(fca)]) is 
proposed, quintessence energy density dominates matter 
energy density asymptotically. 

We conclude that the late time dynamics arising from 
our two-field potential is consistent with the observa- 
tions as well as the theoretical requirements of stability 
of the supersymmetric ground state and the asymptotic 
behaviour of string theory states, provided the observed 
accelerated expansion of the Universe is transient and de- 
celerated expansion is soon resumed, as in Model II. This 
solution has been recently proposed to solve the contra- 
diction between accelerated expansion and string theory 
fjlfl , on general grounds; in this work, we have presented 
a concrete example of a two-field model that exhibits this 
desirable feature. 
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